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Abstract 

The lungs undergo changes that are adaptive for high elevation in certain animal species. In chickens, animals bred at high 
elevations (e.g., Tibet chickens) are better able to hatch and survive under high-altitude conditions. In addition, lowland 
chicken breeds undergo physiological effects and suffer greater mortality when they are exposed to hypoxic conditions 
during embryonic development. Although these physiological effects have been noted, the mechanisms that are 
responsible for hypoxia-induced changes in lung development and function are not known. Here we have examined the 
role of a particular microRNA (miRNA) in the regulation of lung development under hypoxic conditions. When chicks were 
incubated in low oxygen (hypoxia), miR-15a was significantly increased in embryonic lung tissue. The expression level of 
miR-15a in hypoxic Tibet chicken embryos increased and remained relatively high at embryonic day (E)16-20, whereas in 
normal chickens, expression increased and peaked at El 9-20, at which time the cross-current gas exchange system (CCGS) 
is developing. 6c/-2 was a translationally repressed target of miR-15a in these chickens. miR-16, a cluster and family member 
of miR-15a, was detected but did not participate in the posttranscriptional regulation of bcl-2. Around El 9, the hypoxia- 
induced decrease in Bcl-2 protein resulted in apoptosis in the mesenchyme around the migrating tubes, which led to an 
expansion and migration of the tubes that would become the air capillary network and the CCGS. Thus, interfering with 
miR-15a expression in lung tissue may be a novel therapeutic strategy for hypoxia insults and altitude adaptation. 
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Introduction 

The cross-current gas exchange system (CCGS) and blood-gas 
barrier (BGB) morphology are two important factors that 
determine the avian lung air-diffusing capacity [1]. Both play a 
critical role in maintaining embryo homeostasis under low-oxygen 
conditions (hypoxia). Specific processes that take place during 
embryonic lung development and structures in animals that live at 
high altitudes are the basis for highly efficient pulmonary 
ventilation, which leads to altitude adaptation [2,3]. For plain 
chickens at low altitudes, inadequate oxygen exchange results in 
hypoxia syndrome and is lethal [4—6] . Therefore, determining the 
particular altitude adaptation characteristics of high-altitude 
chickens is important. Improving ventilation and protecting 
organs is an important strategy against the high death rate in 
normal animals raised at high altitudes, and studies such as these 
may suggest therapeutic targets for hypoxia syndrome. However, 
the molecular mechanisms of chicken lung development and 
hypoxia adaptation at high altitude have not been well defined. 

For in vitro incubation, chickens are often used to examine the 
interaction between embryonic development and environmental 
conditions. In contrast to mammals, in chicks the mechanisms that 
underlie the development of their lungs involve a complex 
morphology and the generation of a cross-current system. Simple 



lung buds begin to develop at embryonic day (E4), but a complex 
gas-exchanging lung does not develop until E20-21, when the 
chick hatches from the shell. Pulmonary ventilation is an 
important factor for later-stage embryos. Makanya and colleagues 
have found that during chicken embryo development, the atrium, 
infundibulate, and air capillaries (ACs) form in turn from E8 to 
El 9. During development, the ACs branch and anastomose with 
their neighboring cognates during El 6- 19. This process is a direct 
result of epithelial cell attenuation and canalization of the ACs [7]. 
Increasing evidence shows that these processes are accompanied 
by mesenchyme reduction via apoptosis. Lung mesenchyme 
degradation brings the ACs and blood capillaries into close 
proximity with each other [8] . 

Areas of high altitude, such as in Tibet, are regarded as 
"restricted areas of life". The Tibetan chicken, an indigenous 
chicken bred at high altitudes, is well known for its high 
hatchabUity under conditions of reduced oxygen. The hatchability 
of lowland breeds is approximately zero under the hypoxia that is 
simulated at an altitude of ^4000 m, whereas the natural 
hatchability of Tibetan chicks can reach 70-80% [4-6]. An 
oxygen-deprived environment limits embryonic development and 
leads to organ paramorphia [9]. Hypoxia stress adaptation is a 
characteristic of high-altitude animals. But lowland animals also 
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show low-oxygen-induced adaptation during embryonic develop- 
ment. Maina and Kalenga reported that in rat lungs under 
hypoxic conditions, the mean thickness of the BGB is half or less 
than of its thickness under normal conditions [10,11]. Non- 
homogeneous effects of hypoxia in internal organs have been 
observed in avian and mammalian species during development 
[12,13]. Specific weights (organ (lung) weight/body weight) are 
obviously decreased in size from E19 in plain chickens under 
hypoxia, whereas other organs tend to increase in weight [13]. 
Thus, hypoxia induces a lower density lung structure. However, 
the mechanisms of hypoxia-induced embryonic lung development 
and function remain unclear. 

MicroRNAs (miRNAs) are non-coding RNA molecules of 20- 
24 nucleotides long that regulate the expression of other genes by 
inhibiting translation or cleaving complementary target mRNAs 
[14,15]. miRNA functions include the regulation of cell prolifer- 
ation, difiFerentiation, apoptosis, maintenance of stem cell pluri- 
potency, tumor genesis, and others [16-20]. Several miRNAs 
regulate proapoptotic and antiapoptotic genes. miR-15a and miR- 
16, two members of the miR-15a/16 cluster, play a role in 
proapoptosis regulation by inhibiting the translation of the 
antiapoptotic: protein Bcl-2 via binding to the 3 '-untranslated 
region (3'-UTR) of M-2 mRNA [21]. These two miRNAs do not, 
however, always work in the same way [22]. Recent reports have 
focused on the influence of miRNAs on homeostasis and, in 
particular, on the repression of certain genes during hypoxia [23]. 
Because miRNAs have the ability to modify gene expression 
rapidly and reversibly, they are ideal mediators for sensing and 
responding to hypoxic stress. Their ability to modify regulatory 
pathways can affect the ability of an organism to surv'ive under and 
adapt to hypoxic conditions [23]. Therefore, during embryonic 
development, miRNA-mediated regulatory circuits may provide 
flexible and conditional alternatives to embr^'onic development. In 
this study, we identified the influence of long-term hypoxia stress 
on lung development and showed that the miR-15a, downregu- 
lation of HIF-1, was responsive to the oxygen concentration and 
induced mesenchymal ablation through direct inhibition of the 
antiapoptotic gene chicken hd-2 by binding to a unique target 
region. The chick embryo dying may start from an imbalance in 
homeostasis that begins in the lung. The inhibition of miR-15a or 
activation of HIF- 1 or Bcl-2 may prevent hypoxia-induced lung 
damage and reduce chick embryo death. 

Materials and Methods 

Ethics statement 

All animal work was conducted according to the guidelines for 
the care and use of experimental animals established by the 
Ministry of Science and Technology of the People's Republic of 
China (Approval number: 2006-398). The chicken incubation and 
chicken tissue samples collection procedure were approved by the 
Animal Welfare Committee of China Agricultural University 
(Permit Number: XK622). All efforts were made to minimize 
suffering. 

Animals 

The White Leghorn chickens eggs were obtained from the 
Chicken Farm of China Agriculture University. The Tibet chicken 
breed was obtained from Nyingchi Tibet and feeding in the 
Chicken Farm of China Agriculture University. White leghorn 
chicken eggs and Tibet (hT) chickens were divided into four 
groups: hypoxic White Leghorn (hW) chickens and normal White 
Leghorn (nW) chickens, hypoxic Tibet (hT) chickens and normal 
Tibet (nT) chickens. The hT group and hW group were incubated 



under plateau-simulating conditions (37°C; 14% O2; that is the 
similar condition at the 3200 m AMSL (Above Mean Sea Level)). 
The other two groups (nT, nW) were incubated under normal 
condition as control (37°C; 21% O2). Because of the different 
hatchabUity, the starting incubation amount of eggs for each group 
was with big difference. We used 300, .550, 270 and 300 eggs for 
nW, hW, nT and hT samples incubation. Finally, the amount of 
nW, hW, nT, hT samples (eggs/chicks) were 209, 163, 193 and 
192, respectively. The different organ samples (brain, lung, heart 
and liver) were all from these eggs/ chicks. We want to make the 
best value for our egg/ chick samples. 

Northem blot analysis 

Total RNA was isolated from different tissues from adult 
chickens and chicken embryos using TRIzol (Invitrogen, Carlsbad, 
CA, USA). Total RNA (20 |.ig) was fractionated using a 15% 
denaturing polyacrylamide/8M urea gel. The RNA was then 
transferred to a GeneScreen Plus membrane (Perkin-Elmer, 
Waltham, Massachusetts, USA). The 5' ends of the DNA probes 
were labeled with [^^P]ATP (Amersham, Pittsburgh, PA, USA) 
using T4 polynucleotide kinase (New England BioLabs, Ipswich, 
Massachusetts, USA). Hybridization and washing were performed 
as described [24]. Hybridization signals were detected using a 
Phosphor Screen (Molecular Dynamics, Sunnyvale, CA, USA). 

RNA purification and real-time PCR 

Total RNA was isolated with Trizol reagent (Invitrogen) 
ac[:ording to the manufacturer's protocol. After treatment with 
DNase I (RNase-free; TAKARA), oligo(dT) was used to synthesize 
single-stranded cDNA from 2 ^g total RNA using M-MLV 
reverse transcriptase (Invitrogen). miRNA-specific reverse tran- 
scription (RT) stem-loop primers were as follows: 

5sRT: 5'-CTCAACTGGTGTCGTGGAGTCGGCAATT- 
CAGTTGAGAAGCCTAC-3 ' ; 

1 6RT: 5 '-CTC AACTGGTGTCGTGGAGTCGGCAATT- 
CAGTTGAGCACCAATA -3'; 

15aRT: 5'-CTCAACTGGTGTCGTGGAGTCGGCAATT- 
CAGTTGAGACAAACC-3'; 

144R: 5'-CTCAACTGGTGTCGTGGAGTCGGCAATT- 
CAGTTGAGGAGTACA-3 ' . 

miRNA was amplified using SYBR Green (Applied Biosystems, 
Carlsbad, CA, US) and the TaqMan Universal PCR Master Mix 
Kit (Applied Biosystems) with the 7900HT Fast Real-Time PCR 
System (Applied Biosystems). Primers and the TaqMan probe 
were designed using Primer Select software (Primer Express v2.0. 
Applied Biosystems). Real-time primer sequences were as follows: 

GAPDH forward primer: 5'-CGATCTGAACTACATGGTT- 
TACATGTT-3'; 

GAPDH reverse primer: 5'-CCCGTTCTCAGCCTTGACA- 

3'; 

bd-2 forward primer: 5'-AGCGTCAACCGGGAGATGT-3'; 
hcl-2 reverse primer: 5'-GCATCCCATCCTCCGTTGT-3'; 
HIF-1 forward primer: 5'-CAGGTACAA GAGCAACCAA 
CCA-3'; 

HIF-1 reverse primer: 5'-TGGATAAT GACATGGCTA 
ATGAATTC-3' 

HIF-1 probe: 5'-FA]\'I-AGTTCACCTGAGCCC-MGB-3'; 

universal reverse primer: 5'-CTCAACTGGTGTCGTG- 
GAGTC-3' 

miR- 1 44: 5 '-CGCGGCTACAGTATAGATGATG-3 ' 
miR-15a: 5'-GCTGG TAGCAGCACATAATGG-3' 
miR-16: 5'- GCGGAGTAGCAGCACGTAAA-3' 
RN5S: 5'- GCTCTGGAATACCGGGTGCTGT-3' (Table SI 
and S2). 
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Western blot analysis 

Equal amounts of protein were resolved by SDS-PAGE and 
transferred to polyvinylidene difluoride membranes. Membranes 
were incubated with Bcl-2 (C-2) antibody (1:500; Santa Cruz 
Biotechnology, Santa Cruz, CA, US) or HIF-la (C-19) antibody 
(1:500; Santa Cruz Biotechnology). Primary antibody binding was 
visualized with horseradish peroxidase-conjugated secondary 
antibody and detected with enhanced chemUuminescence (Beyo- 
time, Jangsu, China). Secondary antibodies were goat anti— mouse 
IgG (H+L) (1:10000; Santa Cruz Biotechnology) and rabbit anti- 
goat IgG (1:10000; Santa Cruz Biotechnology). For a loading 
control, membranes were reprobed with primary antibody against 
actin (20,1055; 1:5000; Santa Cruz Biotechnology). 

Immunohistochemical staining 

ParafEn-embedded tissues were sectioned (4 [im thick) and 

mounted on glass slides. They were then dewaxed followed by 
antigen retrieval in 0.0 IM sodium citrate at 95' C for 10 min. 
Immunohistochemistry was done using mouse monoclonal anti- 
Bcl-2 (C-2) (IgGl; 1:100; Santa Cruz Biotechnology) in blocking 
reagent at 25°C for 20 min. The primary antibody was omitted for 
negative controls. The slides were then treated consecutively with 
horseradish peroxidase-conjugated goat anti-mouse IgG (1:200; 
Chemicon International, Inc., Teme[:ula, CA, USA) and incubat- 
ed for 1 h at room temperature. Shdes were then incubated using 
a DAB substrate kit (Roche Applied Science, Indianapolis city, IN, 
USA), and the color reaction was allowed to develop for 5-10 min. 
After washing, slides were coversUpped. 

In situ hybridization 

In situ hybridization was performed as described [25]. 
Hybridization was done using the S'-digoxigenin-labeled miR- 
CURY LNA microRNA Detection Probes anti-gga-miR-16 and 

anti-gga-miR-15a (EXIQON, Vedbaek, Denmark). The S'- 
digoxigenin-labeled Scramble-miR Probe (EXIQON) was used 
as a negative control. 

TUNEL Staining 

TUNEL staining was performed using the in situ Colorimetric 
TUNEL Apoptosis Assay Kit (Beyotime), according to the 
manufacturer's manual. Briefly, incubated the frozen sections in 
the cold 0.1% Triton X- 1 00 for 2 min. Then moved the slides into 
the 0.3% H2O2 in methanol for 2 min at RT. Washed 3 times 
with PBS. That was followed by incubation with TdT enzyme 
solution for 60 min at 37°C. The reaction was terminated by 
incubation in a stop/wash buffer for 10 min at 37°C. Then used 
the Streptavidin-HRP for coloration. Incubated the sections in the 
Streptavidin-HRP working solution for 30 min at RT, following 
with the colour reagent incubation for 30 min. 

Cell culture 

The human embryonic kidney 293T cells (American Type 
Culture Collection, Rockville, MD) were cultured at 37°C in 5% 
CO2 in DMEM supplemented with 10% fetal bovine serum 
(HyClone, Logan, UT, US) and antibiotics. 

Plasmid transfection and luciferase assay 

The chicken bcl-2 3'-UTR containing a miR- 15 a— binding site 
and a miR-16-binding site was cloned into the psiCHICK-2 
vector (Promega, Madison, WI, US) using the following primers: 

forward primer: 5'-CTCGAGAGTCACCCAGTTTATCGT- 
3'; 



reverse primer: 5'-CTCGAGGATTCTTCCGCTTCGTCA- 
3'.(Table S3) 

The synthesized fragments containing the reverse complement 
sequences of gga-miR-15a and gga-miR-16 were also individually 
cloned into the psiCHICK-2 vector and used as a positive control. 
The psiCHICK-2 vector was transfected into the human 
embryonic kidney 293T cells with miR- 15a/ 16 mimics (AppUed 
Biosystems) or a miR- 15a/ 16 inhibitor (Applied Biosystems) or a 
negative control (Applied Biosystems). The concentration of mimic 
is 10 pmol. The relative proportion of mimic and inhibitor was 
1:10. The firefly and Renilla luciferase activity of each transfection 
was determined with the Dual-Glo Luciferase Assay System 
(Promega) 44 h post-transfection. 

Results 

The expression of miR-15a and miR-144 is induced by 
hypoxia stress in chicken embryo 

Xu and coUeagues [26] published that the expression levels of 
both miR-15a and miR-144 are higher in the lung tissue. Our 

studies confirmed the result and the tissue special expression. In 
later stages of chick embryo development, relatively strong 
expression of miR- 15a and miR-144 was detected in lung tissues 
with Northern blotting, whereas expression was weak in other 
organs (Fig. lA). To determine whether both miR- 15a and miR- 
144 are highly expressed throughout chicken embryonic develop- 
ment or just a temporally variable expression pattern, we 
compared the miR- 15a and miR-144 expression profiles in lung 
tissues. We set the hatcher to simulate the plateau condition 
(AMSL is 3200 m, 37°C, 14% Og, PO2 is 14.093 kPa) and 
incubated both the Tibet chicken eggs and White Leghorn chicken 
eggs in it to simulate the hT and hW groups. The nT and nW 
chicken groups are the eggs incubated in normal plain condition 
(37''C, 21%02, PO2 is 100 kPa). Lung samples were obtained 
from El 3 to the third day after hatching (d3). Real-time PGR 
showed no difference in miR- 15a and miR-144 expression from 
E13 to E17 in the four groups (data not shown). At E18, the hT 
group showed significant upregulation of miR- 15a that remained 
relatively high and stable through d3 (Fig. IB). miR-15a in hW 
showed continuously increasing expression from E18 and a peak at 
E20 and the hatch out day (dl) (Fig. IB). Expression of miR-15a in 
nW chickens nearly coincided with that of nT chickens except for 
a relative peak value for nW at E19 (Fig. IB). As a control, miR- 
15a in nT chickens showed a relatively smooth line, indicating 
stable, low expression of miR-15a (Fig. IB). A comparison of 
interclass differences showed that at El 8, expression of miR- 15a 
was higher in lung tissues from hT than in the lung tissues from the 
other three groups. At El 9, the time at which the respiratory 
CCGS develops, miR- 15a was also upregulated in hW and nT 
chickens. Comparing the expression levels at E20, we found that 
expression of miR- 15a was relatively low in both tissues under 
normal conditions (nT, nW chickens), and relatively high in the 
hW and hT sample which incubated under plateau simulating 
condition (Fig. IB). Quantitative analysis of miR-144 revealed no 
significant changes over time (data not shown). As a cluster and 
family member, miR- 16 was reported to have the same function as 
miR- 15a [21]. To determine whether the peak value of miR- 15a 
expression was lung specific or was a systemic reaction to hypoxia 
and whether expression of miR- 15a was the same as that of miR- 
16, we performed real-time PGR for miR- 15a and miR- 16 using 
chick embryo lung, heart, brain and liver at E16 and El 9. hW 
showed remarkable upregulation of miR- 15a at E 19 as compared 
with that at E 16 in the embryonic lung (Fig. IC). In the heart and 
brain, miR- 15a expression in hW and hT was also increased from 
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E16 to El 9, but the range was not as large as that in the lung 
(Fig. IC). In the liver, no differences in expression between the two 
embryonic stages or among the experimental groups were 
identified (Fig. IC). We identified a small increase in miR-16 
from E16 to El 9 in embryonic lung in hW (Fig. IC). However, 
there were no clear changes in miR-16 in the other groups. Thus, 
miR-16 also appeared to be affected by hypoxia at a time when 
the respiratory gas exchange was developing. Thus, miR-15a 
expression was sensitive to oxygen concentration, and the strongest 
response was detected in lung tissue. We also found differences in 
hypoxia-induced expression of miR-15a and miR-16 between 
high-altitude and plain chicks. The period of development of 
pulmonary ventilation is crucial for development of lung organ 
function. Although miR-16 belongs to the same cluster and family 
as miR-15a [21], they were differentially sensitive to hypoxia. 
miR-144 also appeared responsive to hypoxia (data not shown), 
but it is not involved in temporal lung development. 

Hypoxia stress upregulates HIF-1a and Bcl-2 protein 

Hypoxia-inducible transcription factor 1 alpha (HIF-la) is 
either anti-apoptotic or pro-apoptotic, according to the cell type 
and experimental design. Because severe or prolonged hypoxia 
induces apoptosis [27], HIF-1 can be used to monitor this 
hypoxia-induced response. To explore the reaction of the lung to 
hypoxia, we detected HIF-la transcripts using quantitative PGR 
and HIF-la protein using western blotting (Fig. 2A, G). Hypoxia- 
induced HIF-la mRNA expression in hW increased from El 6 to 
El 9, whereas HIF-la in hT kept higher (from El 8 to El 9), until 
E20 the trend reverse. In the normoxia experimental groups, nW 
and nT chickens, HIF-la mRNA expression showed a slight 
upward trend, but there was no cUfference between these two 
groups (Fig. 2A). These changes were further confirmed with 
Western blotting, which showed that HIF-la protein was induced 
by hypoxia stress (Fig. 2C). Interestingly, HIF-la protein in nW 
chickens was upregulated compared to nT chickens, especially at 
E19 (Fig. 2C). During normoxia, HIF-la expression can be 
regulated at the translational level [28]. These results suggest that 
nW also suffered acute hypoxia stress at El 9. 

Bcl-2 is downstream of HIF-la in the HIF-la— induced anti- 
apoptotic pathway [29]. Because hypoxic cell death is prevented 
by anti-apoptotic proteins such as Bcl-2 [30], we investigated hcl-2 
mRNA expression during the whole embryo development and 
identified no remarkable variation in bcl-2 mRNA (Fig. 2B, and 
data not shown). However, Bcl-2 protein was strongly decreased in 
E19 hW lung samples as compared with its expression at E16 and 
E18 (Fig. 2C). Moreover, in other groups, Bcl-2 protein levels were 
decreased to different degrees in association with the degree of 
hypoxic stress (the time length) (Fig. 2C). Thus, the regulation of 
Bcl-2 protein levels may not be at the transcriptional level but 
rather at the translational level. We note that miR-15a expression 
was increased at E19 during hypoxia (Fig. IB). 

bcl-2 is a target gene of miR-15a but not miR-16 in 
chicken lung tissues under Fiypoxia 

miR-15a and miR-16 negatively regulate bcl-2 by direcdy 
binding to a particular sequence in the 3'-UTR of bcl-2 and 
inhibiting its translation [21,31]. This repression is sufficient to 
induce apoptosis [32] . Our observation of hypoxia-induced miR- 
15a expression (Fig. IB) and reduced Bcl-2 protein levels at El 9 
(Fig. 2C) indicated an inverse relationship between miR-15a and 
Bcl-2 protein expression, which suggested a causative role for miR- 
15a in the downregulation of hcl-2. However, as a cluster member 
of miR-15a, miR-16 was identified as hypoxia insensitive. To 
further verify whether miR-15a or miR-16 repressed chicken bcl-2 



by binding directiy to the predicted binding site in its 3'-UTR as 
shown in human cells [31], we analyzed the complementation of 
the chicken bcl-2 mRNA (NM_205339.1) and chicken miR- 
15a/ 16 using four different prediction algorithms: TargetScan 
[33], RNA22 [34], DIANA [35,36], and PicTar [37]. Interest- 
ingly, human bcl-2 mRNA (NM_000633) has one target located at 
nt 2529-2536 of its 3'-UTR [21] tiiat is not conserved witii tiie 
chicken bcl-2 3'-UTR. Chicken bcl-2 was predicted to contain a 
miR-15a/16: bcl-2 site located at nt 369-390 of its 3'-UTR 
(Fig. 3A). The RNA22/PicTar results showed that the target site 
sequence complemented only with chicken miR-15a, and the 
TargetScan result suggested that both the chicken miR-15a and 
miR-16 would work (data not shown). We cloned the entire 
chicken bcl-2 3'-UTR, containing the miR-15a/16 binding site, 
and inserted it into the multiple cloning sites of the psiCHICK-2 
vector. We transfected this luciferase reporter vector into 293T 
cells. We also co-transfected these cells with chicken miR-15a 
mimic, mimic and inhibitor, or mimic control (miRNA mimics are 
small, chemically modified double-stranded RNAs that mimic 
endogenous miRNAs and enable miRNA functional analysis by 
up-regulation of miRNA activity. miRNA inhibitors are small, 
chemically modified single-stranded RNA molecules designed to 
specifically bind to and inhibit ('ndogcnous miRNA molecules and 
enable miRNA functional analysis hy down-regulation of miRNA 
activity. The controls are small control molecules). We performed 
the same experiment with chicken miR-16. Chicken miR-15a 
decreased luciferase activity of the reporter vector containing the 
3'-UTR sequence, and the level of luciferase activity recovered in 
the presence of the chicken miR-15a inhibitor. However, no 
change was observed with the chicken miR-16 mimic (Fig. 3C). 
Together, the results suggested that chicken miR-15a decreased 
chicken bcl-2 translation by direcdy acting on a miR-15a-specific 
response element in the 3'-UTR of chicken bcl-2 mRNA and that 
this effect may be different in chicken bcl-2 as compared with that 
in human bcl-2. Chicken miR- 1 6, another member of the cluster, 
did not affect hcl-2 regulation. 

Hypoxia-induced miR-15a overexpression promotes the 
formation of CCGS and a thin BGB 

At E 1 9, the lung begins to function, but the development of the 
complex CCGS is stiU occurring. To further understand a possible 
role for hypoxia-induced miR- 15a in the development of a 
functional lung in vivo, we immunostained sections from hT, hW, 
nT, and nW chicken lungs for Bcl-2 protein and performed 
terminal deoxynucleotidyl transferase nick-end labeling (TUNEL) 
staining. 

Hypoxia upregulated Bcl-2 protein expression at E 1 8 (Fig. 4G, 
H). The strongest staining was of hW chicken lung sections at E18 
(Fig. 4H). Staining of hT and nW lung sections at E18 were 
similar. Their staining were weaker than what detected of the hW 
sections (Fig. 4F, G). No staining was observed in nT chickens at 
E18 (Fig. 4E). As expected, the staining was localized in the 
mesenchyma around the AGs of the chicken lung (Fig. 4F', G', H'). 
No special staining signal was detected in other sections (Fig. 4A, 
B, C, D). 

Apoptotic cells were observed at El 9, except for in nT chickens 
(Fig. 5K). TUNEL staining was darker in hW lung sections at El 9 
(Fig. 5N) and in hT lung sections at E19 (Fig. 5M) as compared 
with staining in nW lung sections at E19 (Fig. 5L). Sections viewed 
at higher magnification showed that apoptotic cells were located in 
the mesenchyme around the migrating epithelial tubes, including 
capillaries (Fig. 5L', M', N'). Interesting, at El 9, no significant 
difference of Bcl-2 protein expression was detected in all the 
groups (Fig. 4K-N). Taken together, these results suggested that 
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Figure 1. Expression of miRNAs in chicken lung tissues. (A) Total RNA from tissues from E19 chicl<en embryos was blotted with probes for mlR- 
15a, miR-144 and U6 (loading control). miR-15a and miR-144 were identified in late-stage embryonic lung tissue. In each sample, the total RNA was 
mixed with samples from 9 chickens. (B) Quantitative expression analysis of miR-15a showed hypoxia-related expression that was affected by both 
the species and environmental conditions. The hW chicken group was most sensitive relative to other three groups. In the nW chicken group, there 
was a response at El 9. At El 9 of the nW chicken group, the expression of miR-1 5a remained relatively high in the hT chicken group compared with 
the nT, nW chicken groups and was largely unchanged in the nT chicken group through the whole embryo stages. Data are expressed as the mean ± 
SEM for each group. (C) Quantitative expression analysis of miR-1 5a and miR-1 6 in the embryonic lung, heart, brain and liver at El 6 and El 9 tissues 
and were expressed as the mean ± SEM for each group. Under hypoxia stress, miR-1 5a was more highly expressed at El 9 than at El 6 in the brain, 
heart and lung for the hW group and in the lung and brain for the hT group. miR-1 6 showed a weak response to stress in the embryonic lung (hW). 
Result statistically different are indicated with an asterisk/s (* P<0.05; ** P<0.01; ns : not significant). E16-20 = embryonic d13-20, respectively. d1, 
d2, d3 = the 1" day after hatching, the 2"'' day after hatching, the 3"' day after hatching. 
doi:1 0.1 371 /journal.pone.0098868.g001 



the lung appeared resistant to apoptosis at El 8, which is likely the 
result of upregulation of Bcl-2. However, at El 9, the time of 
further lung development, the tissues became more sensitive to the 
oxygen concentration, and apoptosis was induced by the hypoxia- 
induced upregulation of miR-1 5a. Of interest, the tube density in 
hW at E19 was also higher than that in sections from those other 
two groups. That indicated adaptive lung development which has 
also been reported in mammals living as altitude adaptation 
[3,38]. These responses to hypoxia stress indicate an oxygen- 



concentration dose effect process in the chick lung. The apoptotic 
cells were located in the mesenchyme around the migrating tubes. 
These mesenchymal cells are in direct contact with the air and are 
the targets of apoptosis, which may lead to expansion and 
migration of tubes, implying precise temporospatial regulation. 
During the subsequent CCGS establishment when the AGs are 
branching and anastomosing with their neighboring cognates, 
mesenchymal cell apoptosis participates in the approximation of 
AGs and blood capillaries (BGs), to establish a thin BGB. 
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Figure 2. HIF-1a and Bcl-2 expression in lung tissue under liypoxia and normoxia conditions. Data were obtained from hi, hW, nT, and 

nW chicl<en tissues at El 6, El 8, El 9, and E20. (A) Quantitative expression analysis of HIF-lot mRNA in embryonic lung at El 6, El 8, El 9, and E20. The 
stress reaction in hW was robust as compared with the smooth response in hT. Data are expressed as the mean ± SEIVl for each group. (B) 
Quantitative expression analysis of bcl-2 mRNA in the embryonic lung, heart, brain and liver at El 6 and El 9. There were no changes in the expression 
of bd-2 mRNA between the El 6 and El 9 in each kind of embryonic tissues. Data are expressed as the mean ± SEM for each group. (C) Analysis of HIF- 
1a and Bcl-2 expression at the protein level in lung tissue (shown in western blot and densitometry value). HIF-la protein expression increased from 
El 6 to El 9, however in different level in hT, hW, nT, and nW group. Bcl-2 protein did show different levels of expression across different time points 
and different groups. (* P<0.05; ** P<0.01; ns = not significant). 
doi:l 0.1 371 /journal.pone.0098868.g002 



Discussion 

Adaptation to hypoxia of the Tibet chicken appears to involve 
numerous changes at the organism's level. These chicken have 
bigger hearts and lungs, higher hemoglobin oxygen transfer 
efhciency and other adaptation characteristics [13,39-42]. Such 
adaptive changes may counteract the reduced oxygen uptake and 
mitigate tissue injury due to hypoxia. There are two oxygen 
sensitive periods during the chick embryo development, the 
embryo day 0-4 (early stage) and the respiratory gas exchange 
(around the El 9, the late stage) [43]. Increased oxygen tension in 
the early and late stages of development results in increased 
hatchabUity, while during while the middle stage of development 
increased oxygen causes little or no change [44^6]. 



The major finding of our study is that, during the late stage of 
chick embryo development, lung mesenchymal cells become 
sensitive to the oxygen concentration when the lung assumes the 
gas-exchange function from the aJlantois at El 9. Hypoxia stress 
induces proapoptotic chicken miR-15a expression and subse- 
quently inhibits the translation of the antiapoptotic protein Bcl-2, 
resulting in mesenchyme ablation. Chicken miR-16, a cluster and 
family member of chicken miR-15a, does not affect the hypoxia- 
induced pathway. In vivo, mesenchymal cell death was seen 
around the tubes, and thus, apoptosis led to expansion and 
migration of new tubes. The mesenchyme ablation led to the 
establishment of CCGS and a thin BGB. Thus, BGB formation 
was associated with the oxygen concentration. The distance 
between AGs and BCs will also be affected by this process. Because 
diflferences exist among species, long-term hypoxia stress in normal 
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Figure 3. bcl-2 is a target gene of miR-15a, but not miR-16 in cKiici^en lung. (A) miR-15a regulates the translation of bd-2 mRNA through a 
sequence that is not conserved with human sequence. miR-15a and the binding site in the gga-bc/-2 3'-UTR are shown, but miR-16 shows no target 
site in this part of the sequence. (B) miR-15a/16 have a consistent target site in human bd-2 3'-UTR. The miR-15a binding site in the bd-2 3'-UTR 
sequence mediates translation repression by miR-15a. (C) The luciferase reporter vector contains two luciferase cDNAs, Renilla luciferase (hRluc) and 
firefly luciferase (hluc). The bd-2 3'-UTR was fused to the hRluc cDNA downstream sequence. In co-transfected cells, the miR-15a mimic decreased the 
expression of hRluc and miR-15a mimic inhibitor rescued hRluc activity; no differences were seen for miR-16. Data are expressed as the mean ± SEIVl. 
* P<0.05, P<0.01 . 
doi:1 0.1 371 /journal.pone.0098868.g003 



chicks is letlial. The lung is the most sensitive organ at later 
incubation times, and an imbalance between organ ablation and 
body development will be lethal (Fig. 6A). To our knowledge, our 
experiments have provided the first evidence that establishment of 
a thin BGB can be mediated by the oxygen concentration during 
the sensitive period of lung development. Development of a thin 
BGB requires precise temporospatial regulation and is tighdy 
correlated with the incubation conditions. The response to stress 
controls mesenchymal cell death via a miRNA-mediated mech- 
anism. We also show that of the two cluster and family members, 



only chicken miR-15a is responsive to hypoxia stress and 
participates in downstream regulation. 

The ability to sense and respond to hypoxia is of fundamental 
importance during embryonic development as well as environ- 
mental adaption. Dysregulation of oxygen homeostasis is an 
important cause of embryo death in hypoxic conditions. An 
organism's oxygen homeostasis begins with the lung gas-exchange 
system. Among the extant vertebrate taxa, the avian lung is 
reported to be the most efficient and complex gas exchanger and is 
a CCGS [47,48]. Microscopy was used to identify the functional 
development of the avian lung in 2006 and showed establishment 
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Figure 4. Immunohistochemical staining (darl< brown coloring) for Bcl-2 protein in lung specimens from chickens incubated under 
normal and hypoxic conditions. Bcl-2 protein expression in lung samples from nW, hT, and hW gradually increased moderately from El 6 to El 8 
before the onset of lung functioning (E-H), and disappeared at El 9. By El 8 in nW lungs, weaker Bcl-2 expression was observed relative to hT (G) and 
hW (H) chicken groups. In nT chicken sections, Bcl-2 protein staining was nearly unchanged throughout the developmental stages examined (A, E, K). 
To see details of the El 8 lung structure, higher-magnification photomicrographs were taken. Bcl-2 staining was detected in the mesenchyme (arrows) 
around the ACs and not in the infundibula or atrias of chicken lung (F', G', H')"- In the hW section at El 8, Bcl-2 staining was strong (arrow in H'), but 
weaker in hT and nW at this stage (arrows in G', F'). Scale bar = 200 |im. 
doi:l 0.1 371 /journal.pone.0098868.g004 



of a thin BGB by the time of hatching [8,49-52]. However, the 
lung of various domestic chicks begins to develop at El 8 [8]. On 
later days (around El 9), the mesenchymal cells of chick lung 
immediately surrounding the migrating epithelial tubes become 
the initial targets of apoptosis. In this study, we report for the first 
time that at El 9, initiation of pulmonary ventilation, which brings 
air into the lung ACs, was induced immediately and promoted 
establishment of an oxygen concentration-dependent CCGS. The 
induced establishment of the CCGS was implemented by induced 
mesenchymal apoptosis, resulting in expansion and migration of 
tubes. The observation that the lung organ/body weight ratio of 
normal chickens was lower than in hypoxic conditions indirectly 
reflects the hypoxia-induced disintegration of mesenchymal cells 
[13]. Adaptive lung development has also been reported in 
mammals living at higher altitudes [38]. E19 is the peak of 
mortahty for embryos growing in both normal and hypoxic 
conditions [53]. Our data showed relative hypoxic stress at El 9 
even when embryos are incubated in normal conditions. We 
suggest that when an acute injury disrupts developmental 
homeostatic regulation, hypoxia may result in the of the affected 
animal. 

miRNAs were discovered nearly two decades ago, and their 
importance was highlighted when they were identified as 
conserved in higher eukaryotes. Recent reports have focused on 



the influence of miRNAs on homeostatic and adaptive develop- 
ment, and in particular, response to hypoxia. An ideal mediator of 
hypoxic stress should be able to modify gene expression both 
rapidly and reversibly, a widely accepted characteristic of miRNAs 
[23]. Therefore, miRNAs may be a molecular rheostat that 
facilitates survival and adaptation to hypoxic conditions by tuning 
and switching regulatory circuits on and off. Various miRNAs 
appear to regulate hypoxia responses in a variety of different 
organisms, cell types, and disease states [22,54—57]. More than 90 
hypoxia-regulated miRNAs (HRMs) have been described, and 
many regulate certain cell types only. However, as a general 
characteristic, these HRMs respond to hypoxic stress rapidly (less 
than 48 h). These hypoxic stress responses are regulated through 
different pathways. miR-107 is directly regulated by p53, which is 
one of the most frequently mutated genes in human cancers [58— 
60]. HIF-1 was identified as target gene of miR-107 both in vitro 
and in vivo, suggesting a p53-miR-107-HIF-l hypoxia response 
pathway [60]. miR-206 was also proved as a direct regulator of 
HIF-1 alpha. In cultured hypoxic pulmonary artery smooth muscle 
cells (PASMCs), miR-206 promotes the Hypoxia-induced pulmo- 
nary hypertension though the HIF-1 a/Fhl-1 pathway [61]. HIF-1 
and microRNAs buUd a network of hypoxia response, HIF-1 was 
reported as regulator of miR-9, it induced upregulation of miR-9 
in PASMCs. As a phenotypic switch, knockdown of miR-9 was 
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Figure 5. Analysis of apoptosis in lung specimens from chickens incubated under normal and hypoxic conditions. From E1 6 to E1 8, no 

TUNEL staining was identified (A-H). At El 9, apoptotic cells were localized in the mesenchyme surrounding the atrias and infundibula of the chicken 
lungs (L, M, N). There was no obvious staining in nT chicken at El 9 (K). At higher magnification, staining was clearly seen in the regions between ACs 
and not in the parabronchi, atrias, or infundibula (arrows in L', M', N'). hW staining at El 9 (arrow in N') was clearly darker than that observed in hT or 
nW lung sections at this stage (arrows in M', N'). The tube density in hW at El 9 was also higher than that in sections from those other two groups. 
Scale bar = 200 |im. 
doi:1 0.1 371/journal.pone.0098868.g005 




Figure 6. Schematic representation of hypoxia-induced gga-miR-15a control of mesenchymal ceil death/ablation during functional 
lung establishment and adaptive development. (A) In chicken lung, hypoxia stress stimulates gga-miR-15a expression, and miR-15a may bind 
to the 3'-UTR of gga-fac/-2 and inhibit its antiapoptotic activity through posttranscriptional gene silencing. Mesenchymal apoptosis then occurs, 
which leads to a functional lung with a thinner BGB or an acute reaction resulting in death in different chicken species. gga-miR-16 shows no 
response to hypoxia, and there is no target site for gga-miR-16 in the 3'-UTR region of gga-bcl-2. (B) Mesenchymal cell death/ablation indeed 
appeared necessary for the formation of a thinner BGB and is promoted by hypoxia (the oxygen concentration). Particular transcriptional factors (TF) 
may respond to the hypoxia stress and induce gga-miR-1 5a expression. For plain chickens, the induced apoptosis would be helpful for the formation 
of a thinner BGB, but it could be lethal if they are exposed to low oxygen concentrations for a long time. (C) However, for high-altitude chickens, 
apoptosis is genetically initiated to form a functional BGB, and it is under control of the whole organism (C). Acute hypoxia: related hypoxia condition 
induced by the CCGS occursing around El 9. Chronic hypoxia: long-term altitude condition with low oxygen concentration. 
doi:1 0.1 371/journal.pone.0098868.g006 
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followed by hypoxia exposure attenuation in PASMCs prolifera- 
tion [62]. miR-210 and miR-373 were discovered in HeLa and 
MCF-7 cells. Both miRNAs contain hypoxia-responsive elements 
(HREs) in their promoter regions, suggesting direct regulation by 
HIF-1 [63]. Most hypoxic injury results in cell apoptosis. 
Hayashita and Frankel reported that miR-21 and the miR- 17-92 
cluster function as antiapoptotic miRNAs in various cancers 
[64,65]. In contrast, Cimmino et al. showed that miR- 15a and 
miR- 1 6 promote apoptosis by posttranscriptional gene silencing of 
bcl-2 [21]. In this study, our data suggested that the expression of 
miR- 15a, but not miR- 16, is sensitive to the oxygen concentration 
and was significantly increased in lung mesenchymal cells in 
chicken. This results in select apoptosis in mesenchymal cells 
around the migrating tubes through reduction of Bcl-2 levels, as 
shown in vitro and vivo. Although miR- 15a and miR- 16 belong to 
the same cluster and miRNA family, Yin et al. also reported that 
miR- 15a plays a causative role in the regulation of apoptosis by 
directly targeting bcl-2 during ischemic vascular injury [22]. 
Together, our experimental data extend these early findings and 
further show that chicken miR- 15a can repress bcl-2 translation by 
directly binding to the 3'-UTR of hcl-2 with a different target site 
sequence. To our knowledge, our work is the first to document 
that chicken miR- 15a, but not miR- 16, regulates apoptosis by 
direcdy targeting chicken bcl-2 at a specific target region during 
later chick lung CCGS development. 

Kang et al. [66] reported that the miR-15a/16-l cluster is a 
homo-cluster. Because the homologous miRNA in the homo- 
cluster has almost the same proximal targets or genes in the same 
pathway, homo-cluster miRNAs show a direct regulatory coordi- 
nation involving one step. Thus, there is no intermediate regulator 
between the miRNAs and their targets; direct coordinated 
regulation has the advantage of accuracy and quickness. Hua 
and colleagues identified a group of regulator)' miRNAs, including 
miR- 16, which regulates the expression of vascular endothelial 
growth factor [67]. miR- 16 may function during BC network 
formation, which is perpendicular to the long axis of AGs. miR-16 
shows unequal expression among tissues [68]. In this study, we 
suppose that in chick, members of the miR-15a/16-l cluster 
respond to hypoxia through dilferent target genes and pathways. 
We showed relatively high expression of miR- 144 in chicken 
embryo lungs. Analysis of miR- 144 showed a lower response to 
hypoxia. miR- 144 expression has been reported to be limited to 
the blood islands at stage 8 [69], and it is crucial for erythroid 
homeostasis [70]. miR-144 also modulates oxidative stress 
tolerance and is associated with anemia severity in sickle cell 
disease [71]. Thus, we suggest that miR-144 may be responsive to 
hypoxia but may have no direct relationship with lung develop- 
ment. 

Lung organogenesis is a field that has attracted many 
investigators. Since "molecular embrs'ology" ^vas reviewed, more 
challenges have emerged in developmental lung biology. miRNAs 
play an important role not only in cell proliferation and 
tumorigenesis but also in apoptosis that is associated with organism 
and organ development. Yang et al. identified 167 differentially 
expressed miRNAs during rat lung organogenesis, and Dong 
reported 117 miRNAs that are dynamically and temporally 
regulated during mouse lung organogenesis. Mujahid et al. found 
that, during the important window in lung development, 25 
miRNAs showed different expression profiles across gestation, and 
37 miRNAs changed between males and females [72-74]. 
Transcription factors connect the extracellular environment to 
the signaling network in cells. HIF-1 is a widely recognized 
hypoxic stress transcription factor [23]. In a recent report, HIF-1 
was shown to function in lung branching morphogenesis and 



epithelial maturation. Functional loss of HIF-1 results in neonatal 
lethality and severe respiratory distress. Interestingly, some 
abnormalities were detected in the endothelial barrier, which 

indicates a disorder in development of the epithelial and 
mesenchyme [75]. By controlling both O2 delivery and utilization, 
HIF-1 functions as a important regulator of O2 homeostasis. 
Under tissue hypoxia or ischemia, HIF-1 protect the heart against 
injury by inducing angiogenic growth factors, which stimulate 
vascular remodeling to increase blood flow following [76]. 
Intrestingly, for the High altitude adaptive species, during the 
Tibet chicken embryo lung development, HIF- 1 and other factors 
induce complex CCGS and a thin BGB which is necessary for the 
High altitude adaptive lung fianction. In the present study, we 
detected upregulation of both HIF-1 and miR- 15a with a similar 
expression profile during E16-E19. We hypothesize that the 
induced ablation of mesenchymal cells may involve HIF-1 as a 
hypoxia sensor, which then indirecdy decreases Bcl-2 protein 
levels by inducing miR- 15a expression. However, more studies are 
needed to elucidate the details. Other transcription factors are 
likely to participate in this regulatory network. The sensor of 
oxygen concentration is important here, because apoptosis occurs 
only around the ACs and not at the top to the atria. In the Tibet 
chicken, induced apoptosis is necessary for establishment of 
functional lungs containing complex CCGS and a thin BGB. 
However, when the induced apoptosis is uncontrolled and organ 
homeostasis is lost, normal chickens ^\ ill die for pressure overload 
(Fig. 6B, C). Other molecules have been reported as miR- 15a 
regulators in diflerent cell types or tissues, and some are hypoxia 
regulated. miR-709 responds to apoptosis stimuli and prevents 
expression of miR- 15a/ 16 in mouse [77]. Myc and HDAC3 are 
regarded as miR- 15a/ 16 regulators in mantle cells and other non- 
Hodgkin B-cell lymphomas [78,79]. HERA16, a clinically 
important oncogenic isoform of HER2, is another regulator of 
miR-15a/16 in breast tumors [80]. In human tumors, E2F1 is a 
positive regulator of miR- 15a and miR-16, but only miR- 15a 
inhibits expression of cyclin E [81]. In a mouse model of ischemia- 
induced cerebral injury, PPAR5 (peroxisome proliferator-activated 
receptorS) was identified as a regulator of miR- 1 5a that inhibits 
the induction of apoptosis by preventing miR- 15a expression [22]. 
Intrestingly, three isoforms of PPARs (Qt,Yand 5) were all 
documented direct downstream of miRs [22,82,83]. PPRAo! were 
identified had specific functional variants between highland 
Tibetan populations and normal [85] . However, although PPRAa 
is in the Gene Ontology biological process 'response to hypoxia', 
but it is not even in the HIF pathway [84] . So, additional studies 
win be needed to further understand the interaction between 
miRNAs, gene functional variation and conditional development. 

In this study, we identified the influence of long-term hypoxia 
stress on lung development. In conclusion, the mechanisms of 
hypoxia-induced, miR- 1 5a-mediated CCGS and thin BGB 
estabhshment are important for understanding adaptational 
development of chick lung at particular stages. HIF-1 may 
regulate miR- 15a in this setting. We also showed that only miR- 
15a and not miR-16 was responsive to the oxygen concentration 
and induced mesenchymal ablation through direct inhibition of 
the antiapoptotic gene chicken bcl-2 by binding to a unique target 
region. The stage at which embryos die may results from an 
imbalance in homeostasis that begins in the lung. Therefore, either 
inhibition of miR- 15a or activation of HIF-1 or Bcl-2 may 
potentially be therapeutic options for preventing hypoxia-induced 
lung damage and may reduce death at particular stages of chick 
embryo development. 
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